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Abstract 


Recently, various alternatives to batteries, such as microfabricated fuel cell systems, have been proposed for portable power generation. In 
large-scale power production plants emphasis is placed on energy conversion efficiency. On the other hand, the intrinsic design objective for 
portable power generation devices is the energy density, i.e., the electrical energy generated from a given mass or volume of device and fuel 
cartridge. It is plausible to stipulate that an increase in the energy conversion efficiency of a system leads to an increase in energy density, but we 
demonstrate through theoretical analysis and case studies that the two metrics are not equivalent. In some cases, such as systems with a combination 
of fuels, maximizing efficiency leads to drastically different design, operation and performance than maximizing energy density. Another interesting 
observation is that, due to interaction between components, maximal component efficiency does not always imply maximal system efficiency. 


© 2006 Elsevier B.V. All rights reserved. 
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1. Introduction 


The increasing need for man-portable power generation for 
military and civilian applications is well established, e.g., [1,2]. 
The predominant technology is currently batteries, but there are 
technological, economical and ecological concerns. One of the 
problems is that the energy density, i.e., the energy produced per 
mass or volume of the battery, is relatively low, in the order of 
a few hundred Wh1~! and Wh kg™! for rechargeable batteries 
[3,4]. 

One of the promising alternative technologies is the use of 
microfabricated fuel cell-based systems with common fuels or 
chemicals, such as hydrocarbons or alcohols, as the energy 
source. The promise of this approach is that the energy density 
of these fuels is relatively high and fuel cells can, in principle, 
achieve high efficiencies [5,6]. Moreover, similarly to batter- 
ies, fuel cell systems have few or no moving parts, and run 
silently. In order to achieve portability the use of microfabri- 
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cated devices, as opposed to conventional devices, is plausible 
[7]: 

The area of man-portable power generation is extremely 
active and there are several academic and commercial programs 
exploring microfabricated fuel cell systems. The vast majority of 
the literature deals with fabrication issues, which are outside the 
scope of this paper. There are a few contributions on scaling and 
system-level considerations [8-11], and some contributions on 
detailed modeling [12,13]. The reader is referred to the review 
articles by Holladay et al. [14] and by Maynard and Meyers [15], 
as well as the collection [16]. 

The focus of this paper is a thorough discussion of the 
appropriate design objective for portable power generation. We 
demonstrate that the design objective of maximal energy den- 
sity is not necessarily equivalent to the objective of maximal 
energy efficiency, in the sense that the two objectives may lead 
to drastically different design and/or operation. We first analyze 
several cases for which the two objectives are not equivalent from 
a mathematical point of view and then we demonstrate through 
case studies that this difference indeed affects the optimal design 
and operation. 
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Nomenclature 

Cpi molar heat capacity of species i (J mol~! K7!) 

grav gravimetric energy density (Wh kg7!) 

vol volumetric energy density (Wh17!) 

F Faraday constant (C mol`!) 

Gi? molar gas phase Gibbs free energy of pure species 
i at reference pressure (J mol`!) 

H? molar enthalpy of species i, ideal gas (J mol™!) 

HË molar enthalpy of formation of species i (J mol~!) 

Hh heating value of species i (J mol!) 

AH"? vaporization enthalpy of species i (J mol~!) 


A,H enthalpy of reaction r (J mol`!) 

I current (A) 

J set of species 

J set of elements 

M mass (kg) 

MV; molar volume of species i (1 mol`!) 
MW; molecular weight of species i (kg mol~!) 
Ni, j molar flowrate of species i of stream j (mol s~!) 
P pressure (bar) 

Pref reference pressure (bar) 

PW power (W) 

R gas constant (J mol`! K7!) 

T temperature (K) 


Tamb ambient Temperature (298 K) 

Tref reference Temperature (298 K) 

U voltage (V) 

Uloss overall heat transfer coefficient (W m~? K7!) 
V volume (m?) 

Z number of electrons exchanged 


Greek letters 


di, j number of atoms of element j in species i 

€ product of emissivity and view factor 
(W m?/W m?) 

E tolerance 

n conversion efficiency 

nsorc SOFC efficiency (W/W) 

C conversion of reaction r (mol s~!/mol s~!) 

[hj KKT multiplier corresponding to element j 
(Jmol!) 

Vri stoichiometric coefficient of species i in reaction 
r 

&, extent of reaction r (mol s~!) 

pi density of species i (kg m7?) 

OA maximal allowable tensile stress (Pa) 

T residence time (s) 

Tmission Mission duration (time between refueling) (h) 

@ air ratio (mol s~!/mol s~!) 


2. Analysis of design objectives 


In large-scale power production, emphasis is placed on effi- 
cient utilization of the fuel. This is because the fuel cost is of the 


same or higher order of magnitude as the fabrication cost of the 
power production system. In man-portable power production the 
economical and ecological operating costs are much smaller rel- 
ative to the fabrication costs of the systems. Typically, different 
man-portable power generation systems are compared using the 
metric of energy density of the system [6]. The specific energy, 
or gravimetric energy density Corey [Wh kg7!], is expressed as 
the electrical energy produced per unit mass of system [3] and 
the (volumetric) energy density Bea [Wh1~'] is defined as the 
electrical energy produced per unit volume of the system 
TmissionP W sys TmissionP W 


Evol = vsys á 


sys 


ĉ grav = Msys 


where the mission duration Tmission [h] is the time between refuel- 
ing or recharging, PW [W] is the power output (assumed constant 
for simplicity), M®S [kg] is the mass of the system, and V53 
[1] is the volume of the system. Depending on the application, 
either of the densities is more important. It is essential to define 
the system appropriately including the power generation devices 
as well as the fuel containers. The most promising application 
of fuel cell systems is for long mission durations [9,17]. For the 
limit of an infinite mission duration the system energy density 
approaches the fuel energy density, i.e., the power produced per 
unit massflow or volumetric flow of the fuel 


eluel — PW etuel = PW 
sev 36005; MW;N;in =“ 36009; MV: Nain’ 


where N; in [mol s—!] is the inlet molar flowrate of species i, 
MW; [kg mol—!] is the molecular weight of species i, MV; 
[I1mol~!] is the molar volume of species i at storage conditions, 
3600 is the conversion factor from hours to seconds, and the 
summation is taken over all stored fuels. 

Power generation is associated with heat generation, 
inversely proportional to the overall system efficiency, e.g., [18]; 
inefficient processes might be considered uncomfortable for 
portable applications because of the large heat generation, e.g., 
acellular phone getting hot, or yield an undesired heat signature 
in the battlefield. 

The overall energy conversion efficiency n°% for a power 
system is given by the power generated divided by the chemical 
energy fed to the system per unit time, accounting for unreacted 
fuel 
n = e (1) 

D i H; i N i,in 
where Hheat [J mol~'] is the heating value of species i, i.e., the 
heat released from complete combustion of one mole of species 
i at reference conditions Tief = 298 K. 

Often it is important to define the energy conversion effi- 
ciency of components of a system, e.g., the reforming reactor or 
the fuel cell, and sometimes alternative definitions of efficiency 
are plausible. For instance, the efficiency in the fuel cell can be 
defined as the quotient of power produced to chemical energy 
fed to the fuel cell per unit time or as the quotient of power 
produced to chemical energy consumed per unit time. One defi- 
nition accounts for unreacted fuel while the other does not. The 
latter definition is more appropriate when the fuel cell effluents 
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are utilized in a burner. As is shown in the following for the case 
of one chemical reaction the fuel cell voltage is an appropriate 
metric for the energy conversion efficiency. The definition of the 
component efficiency gives 


_ PW 

E ArH ( Tamb) : 
where nrc is the fuel cell efficiency, PW [W] is the power, & 
[mol s~!] is the extent of the reaction, and A, H(Tamp) [J mol~!] 
is the corresponding enthalpy of reaction at ambient conditions. 


On the other hand, the power produced is given by the product 
of voltage U [V] and current J [A] 


PW = UI. 


NFC 


The current is the product of the extent of reaction, the Faraday 
constant F [C mol" !], and the number of electrons z exchanged 
in the oxidation of hydrogen (2 per molecule of H2) 


I = zF, 


and therefore 
Fz 
A,H (Tamb) i 

The system energy conversion efficiency is a function of the 
component efficiencies. For components in series the overall 
efficiency is the product of the individual efficiencies; for com- 
ponents in parallel the overall efficiency is the average of the 
component efficiencies weighted with the individual power lev- 
els. As is demonstrated in a following case study, improving 
an individual component efficiency does not necessarily imply 
that the energy conversion efficiency of the system will also 
be improved. This is due to the often complicated interaction 
between individual components. 

The objective of maximal energy density is in general not 
equivalent to the objective of maximal efficiency. The simplest 
example illustrating this is a comparison between different fuels; 
choosing a fuel with high energy density can lead to a higher 
system energy density despite a lower efficiency. For instance 
a 35% efficient butane system has a higher energy density than 
a 70% efficient ammonia system (see Fig. 1 in [7]). A similar 
behavior is seen for systems with a combination of different 
fuels/chemicals, where energy density and energy conversion 
efficiency bear different weights on each species. Note that a 
combination of fuels leads to performance improvements in 
some cases [17]. The extreme case of species combination is 
the addition of water in steam reforming reactions, which does 
not (directly) affect the energy efficiency but greatly reduces the 
energy density, assuming that the water cannot be recycled but 
has to be carried along with the fuel. Another extreme example 
is the case that no ambient air is available for the oxidation reac- 
tions, e.g., for underwater or space applications, and an oxygen 
or air cartridge must be carried. The gases require significant 
storage volumes and cartridge mass, and the energy density is 
significantly lower compared to the use of ambient air [7]. On 
the other hand, the energy stored in the compressed gas is very 
small compared to the chemical energy of the fuels and therefore 
the metric of energy conversion efficiency does not account for 
the stored gas. Moreover, the use of pure oxygen may increase 


nsorc = U 


the energy conversion efficiency compared to the use of ambi- 
ent air. For systems involving only one stored species i, the fuel 
energy density and efficiency are proportional, since 


eluel = _ Pw eluel = iy 
grav — 3600MW;, Niin > vel 3600 X; MV; N; in 
po 
A Ni, ih 
and therefore 
f , 
a a E 


nS 3600MW; ° nS 3600MV; ` 


On the other hand, the system energy density is also a function 
of the device size, while the energy conversion efficiency is not. 


3. Case studies 


In this section we use different case studies to demonstrate 
that the design objectives of maximal energy density and maxi- 
mal energy efficiency are not only different from a mathematical 
point of view, but may also lead to different optimal design 
and/or operation and therefore different performance. We first 
use the methodology for comparison of alternatives developed 
in [7,17], which is based on mass and energy balances along with 
prespecified parameters describing component performance. 
Then we use a kinetic-based case study based on the methodol- 
ogy for optimal design and operation of fixed process structures 
developed in [19]. Finally, we present an equilibrium-based case 
study. The three case studies are in a sense different extreme 
cases: the system-level case study allows the comparison of dif- 
ferent fuels and different processes, but has some parameters that 
are not calculated from first principles; the results of the kinetic- 
based case study are limited by the effect of temperature on 
reaction rates and heat losses and are therefore scale-dependent; 
the results of the equilibrium-based case study are independent 
of chemical kinetics and heat losses. The reason for the vari- 
ety of case studies is to demonstrate that the results are not an 
artifact of the methodology, but rather an intrinsic difference 
of the design objectives. Moreover, using different formulations 
we investigate the suitability of the design metrics for the selec- 
tion of a process among alternatives (first case study) as well as 
for the optimal design operation of a given process (second and 
third case studies). 


3.1. System-level case studies 


In [7,17] a methodology for the comparison of micropower 
generation alternatives is presented. It is based on relatively sim- 
ple models that are general enough to be independent of design 
details, such as the choice of catalysts and reactor configura- 
tion. Mass and energy balances are combined with user-specified 
parameters that describe the performance of the system compo- 
nents. Since heat integration and heat losses to the ambient are a 
crucial part of the design, they are considered in parallel with the 
flowsheeting options. This is done by combining units in con- 
ceptual stacks of uniform temperature and by calculating the 
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heat losses to the ambient based on the device surface area and 
constant conductive and radiative heat transfer coefficients. The 
energy balance is closed by oxidizing sufficient fuel in the burner 
units. Out of the thousands of possible process configurations 
and layouts in [17], here we compare four alternatives in terms 
of the achievable energy densities and energy efficiency. For 
the model details, including the calculation of fuel and system 
energy density the reader is referred to [7,17]. 

For all four alternatives we consider a power generation 
of PW = 10W. For the calculation of system energy density 
a mission duration of Tmission = 30h is assumed. Heat losses 
are accounted for with an overall transfer coefficient Uloss = 
3 W m~? K! and an emissivity (including the view factor) of 
€ = 0.2. The operating pressure is assumed to be | atm for all 
four processes. Three out of the four processes employ atmo- 
spheric air. For these processes a pressure increase mechanism 
such as a microblower is assumed, with a power requirement 
of 3 kJ mol`! of air fed to the system. This conservative power 
requirement is calculated assuming isothermal compression at 
ambient temperature from 1 to 1.2 atm with an efficiency of 
15% [7]. The SOFC and burner operating temperatures are set 
to 1000 K and the outlet temperatures to 600 K. The conversion 
is assumed to be 90% in the reactor, 80% in the fuel cell and 
95% in the burner. The residence time in the SOFC is assumed 
to be 20 ms and in the burner 1 ms. The effect of packaging and 
thermal insulation are accounted for by calculating the device 
volume as 10 times the inner volume of the reactors and fuel cell. 
A device density of 1 kg1~! is assumed. The fuel cell power out- 
put is assumed to be 70% of the product of extent of reaction 
and standard Gibbs free energy of reaction calculated at the fuel 
cell temperature 


PW =0.7 5 = 5 viG;°(T), 


where &, is the extent of electrochemical reaction r, v, is the sto- 
ichiometric coefficient of species i in reaction r, and G;°(T) is 
the molar gas phase Gibbs free energy of pure species i at the ref- 
erence pressure. This simplification eliminates the dependence 
on composition. 

The first alternative considered is the combination of ammo- 
nia cracking for the generation of hydrogen with a polymer 
electrolyte membrane (PEM) fuel cell (see Fig. 1). While ammo- 
nia is very toxic, it is considered as a hydrogen source, e.g., [20], 
because the cracking products do not contain carbon monox- 
ide which is a poison for PEM fuel cells. To overcome heat 
losses and the endothermicity of ammonia decomposition, the 
fuel cell effluent from the anode compartment is oxidized. Note 
that using a new air stream for the hydrogen burner leads to a 
better performance than using the air excess from the fuel cell 
effluent of the cathode compartment, because of the low tem- 
perature of the PEM and the nitrogen and water that are present 
in it. The hydrogen burner and the cracking reactor are consid- 
ered to be in an isothermal stack for better performance [17]. 
The residence time in the reactor is assumed to be 1 ms, the 
operating temperature 1000 K and the outlet temperature 600 K, 
accounting for partial heat recovery. Atmospheric air is used for 
the oxidation reactions with an energetic penalty as described 


aste 
PEM 


Fig. 1. Flowsheet for ammonia cracking. 


above. The equilibrium concentration of ammonia at 1000 K is 
approximately 200 ppm and therefore a significant performance 
decrease of PEM may be observed [21,22]. Our models do not 
account for this degradation but for actual deployment this issue 
has to be addressed. A possibility is to use an ammonia sor- 
bent, e.g., [23]. Including the weight and volume of this sorbent 
would affect the system energy density but not the conversion 
efficiency, thus strengthening our point that energy conversion 
efficiency is not the appropriate metric. We also neglect any need 
for cooling or humidifying of the PEM. 

The second alternative considered is methane oxidation in a 
direct solid-oxide fuel cell (SOFC) using compressed oxygen 
for the oxidation (see Fig. 2). At the macroscale steam reform- 
ing of natural gas is the predominant technology for hydrogen 
production. Methane as a feedstock for hydrogen production has 
the advantage that no carbon bonds need to be broken and the 
ratio of hydrogen-to-carbon atoms is maximal among the hydro- 
carbons. Moreover, compared with methanol and formic acid, 
methane has the advantage that it is not already partially oxi- 
dized. As a feedstock for portable applications it has the major 
drawback that it is supercritical at ambient temperatures and 
most likely to be stored as a compressed gas, which leads to 
low energy densities. The ideal case of direct electrochemical 
oxidation of methane in a SOFC is considered. As a further 
idealization we ignore the effect of carbon deposition on the 
catalysts. We assume gas storage in a plastic container with a 
maximal stress Omax = 100 MPa, a density of po = 1.5kg m~? 
and a storage pressure of 10 MPa. 

The third alternative considered is steam reforming of 
propane and oxidation of the generated hydrogen and carbon 
monoxide in a SOFC (see Fig. 3). Propane has the advantage 
that it has a high intrinsic energy density and can be stored as a 
liquid under moderate pressure. Moreover, the vapor pressure of 
propane is sufficient to overcome pressure losses and no com- 
pressor or pump is needed for the fuel. Atmospheric air is used 
for the oxidation reactions with an energetic penalty as described 
above. To overcome heat losses the fuel cell effluent is oxidized 
and since the SOFC operates at a high temperature the air excess 


Fig. 2. Flowsheet for methane oxidation in a SOFC. 
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waste 


waste 


waste 


Fig. 4. Flowsheet for partial oxidation of propane. 


from the fuel cell is used for the burner. The residence time in the 
reactor is assumed to be 2 ms, the operating temperature 1000 K 
and the outlet temperature 1000 K. A stoichiometric mixture of 
propane and water is used. The power generated in the SOFC 
is based on equivalent hydrogen production [17,24], i.e., it is 
assumed that CO and C3Hg undergo internal reforming and that 
only hydrogen is electrochemically oxidized. 

The fourth alternative considered is partial oxidation of 
propane and oxidation of the generated hydrogen and carbon 
monoxide in a SOFC (see Fig. 4). Atmospheric air is used for 
the oxidation reactions with an energetic penalty as described 
above. The exothermicity of the partial oxidation reaction suf- 
fices to overcome the heat losses. The residence time in the 
reactor is assumed to be | ms, the operating temperature 1000 K 
and the outlet temperature 1000 K. The power generated in the 
SOFC is again based on equivalent hydrogen production [17,24]. 

The results of the comparison are summarized in Table 1. 
There are several interesting points to note. The partial oxidation 
of propane leads to the highest energy densities, despite the 
lowest energy conversion efficiency, while the direct oxidation 
of methane leads to the highest energy conversion efficiency 
but the lowest system energy densities. Both the system and 
fuel volumetric energy densities of methane direct oxidation are 
very low, due to the large volume required for the storage of 
the gases (methane and oxygen). Moreover, the gravimetric fuel 
energy density of methane is low, due to the mass of oxygen 
which is accounted for [17]. The large difference in gravimetric 
energy density of methane direct oxidation is due to the mass of 
the gas cartridges. Despite lower energy conversion efficiency, 


the propane-based processes lead to higher energy densities than 
the ammonia-based process, because of the intrinsic difference 
in energy density between the two fuels. The comparison of the 
two idealized possibilities of propane fuel processing reactions 
shows that the energy conversion efficiency is not a suitable 
metric for man-portable applications, because it does not account 
for the water weight and volume; the higher energy efficiency of 
steam reforming is due to the generation of additional hydrogen 
in the reactor. With the exception of the process requiring gas 
storage, the fuel energy density and system energy density give 
the same qualitative comparison among processes. 


3.2. Kinetic-based study 


In this case study we consider the ammonia-butane process 
from [19]. The power generation process consists of a fuel pro- 
cessing reactor, a solid-oxide fuel cell (SOFC) and two burners, 
fabricated in a single silicon stack fed with ammonia and butane 
fuels (Fig. 5). Ammonia is first catalytically decomposed into 
nitrogen and hydrogen. The produced gases are fed into the 
anode of the SOFC. An air stream is fed to the cathode of the 
SOFC. The anode and cathode effluents are finally mixed and 
fed into burner I, along with potentially a second air stream, 
for catalytic oxidation. In parallel, a mixture of butane (C4Hj0) 
and air is fed into burner II for catalytic oxidation to produce 
heat, thus maintaining the stack at a desired, sufficiently high 
temperature, despite the heat losses and the fact that ammonia 
decomposition is an endothermic reaction. 

The energy efficiency is given by the power output divided 
by the product of heating value and molar flowrates (cf. Eq. (1)). 
Based on the heating values and molecular weights of butane and 
ammonia, maximizing energy conversion efficiency for a given 
power output is equivalent to minimizing Nyy; + 8.4Nc,Hyo 
whereas maximizing the energy density is equivalent to mini- 
mizing Nn; + 3.4Nc,4H,,- Indeed, 

fuel 


PW 
“ery 3600 55, MW; Niin 
PW 
3600(0.017 NNH, + 0.058Ncyti;9) 


and since PW is a constant 


max el <> min(0.017NNH, + 0.058 NcyH,9) 


grav 
> min(Nnu, + 3.4Nc4Hyo)- 


The heating values of butane and ammonia can be calculated to 
approximately 2.7 x 10° and 3.2 x 10°Jmol~!, respectively, 


Table 1 

Results for the comparison of processes 

Performance metric NH3 CH4 C3Hs ref. C3Hg POX 
Energy efficiency (%) 27 30 23 16 
Volumetric fuel energy density ( Wh17!) 1570 170 1740 1920 
Gravimetric fuel energy density (Whkg~') 2580 1610 2540 3900 
Volumetric system energy density (Wh 17!) 1510 130 1650 1800 
Gravimetric system energy density (Whkg~!) 2370 320 2320 3340 
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Fig. 5. Conceptual process flowsheet. 


by the oxidation reactions and therefore 
_ PW _ PW 
SD, APN; in 3.2 x 105Nyp, + 2.7 x 106 Nc,Hio 


sys 
ney 


Since PW is a constant 
max nS < min(3.2 x 10° Nnn; + 2.7 x 10°NcyHyo) 
< min(Nyn; + 8.4NcyHio)- 


It is clear that the two objectives bear a different weight on 
each fuel flowrate and running the optimization problems shows 
that the respective optimal designs and operations are drastically 
different. 

We use the same model as in [19], wherein the ammonia 
decomposition reactor, the solid-oxide fuel cell, and the fuel 
cell residual burner are modeled as isothermal and isobaric 
plug-flow reactors. The gas phase is assumed ideal, which is 
plausible because of the low pressure and high temperature. 
The change in density due to reaction is accounted for. For 
the kinetics of ammonia decomposition the reduced one-step 
expression from [25] is used. For the SOFC the kinetic data are 
taken by Achenbach [26]. For the oxidation of fuel cell efflu- 
ents, the kinetic mechanism and data proposed by Pignet and 
Schmidt [27,28] are used. Nitrous oxide (N20) production is 
not accounted for and the unimolecular decomposition of NO is 
neglected. A lumped model, based on global mass and species 
balances, is formulated to describe butane catalytic combustion 
in burner II. It is assumed that the combustion reaction takes 
place to a fixed conversion and an excess of oxygen is always 
considered. 

The optimal design and operation of the system is formulated 
as a nonlinear optimization problem with differential-algebraic 
equations embedded. The optimization variables are the sizes of 
the ammonia decomposition reactor, fuel cell and burner (design 
variables) as well as the fuel flowrates, air flowrates and fuel cell 
voltage. A common energy balance is considered for the device, 
accounting for heat losses to the ambient. Note that the fuel cell 
voltage affects both the power production and heat generation. 
We consider a power output of PW = 10 W with an operating 
temperature T = 1000K, an outlet temperature Tout = 650 K 


and the same values for the various design parameters as in the 
base case of [19]. 

Table 2 summarizes the key findings for the two objective 
functions. When maximizing efficiency the achievable effi- 
ciency is approximately 25% higher compared to the case of 
maximizing energy density at the expense of a 10% reduction 
of the energy density. As expected from the weighting of the 
two fuel flowrates, maximizing the energy efficiency leads to a 
smaller butane flowrate at the expense of the ammonia flowrate. 

A very interesting result is that the optimal device size is 
significantly smaller in the case of optimal efficiency, mostly 
due to a decrease in the fuel cell volume. Note that in the case 
of maximal energy density butane is not penalized as much and 
therefore the optimal design allows for a bigger device with a 
better hydrogen conversion in the fuel cell but higher heat losses, 
resulting in higher butane flowrates. On the other hand, in the 
case of maximal energy efficiency the optimal design is such 
that the ammonia line is nearly autothermal and a small butane 
flowrate is required to close the energy balance. 

The fuel cell voltage in the case of maximal efficiency is sig- 
nificantly lower than in the case of maximal energy density. The 
most likely explanation is that in the case of optimal efficiency 
more heat generation is required from the fuel cell, which results 
in a lower efficiency and therefore lower voltage. This result also 
emphasizes the point that maximization of the efficiency of one 
component does not necessarily lead to maximization of the 
system efficiency. 


Table 2 
Performance parameters for kinetic-based case study 


Property Maximal energy Maximal energy 
density efficiency 

Energy density ( Wh kg7!) 1247 1107 

Energy efficiency (%) 17 21 

Fuel cell conversion (%) 92 74 

Fuel cell voltage ( V) 0.43 0.33 

Fuel cell efficiency (%) 44 33 

Device inner volume (cm?) 10.5 6.9 

Ammonia flowrate ( mol s7!) 89 145 

Butane flowrate ( wmol s7!) 12 0.8 
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3.3. Equilibrium-based case study 


In the last case study we consider equilibrium-limited pro- 
duction of hydrogen from octane via steam reforming and partial 
oxidation. We neglect heat losses, assuming that the power level 
is sufficiently high to justify this. This case study is therefore 
mostly relevant for applications in the order of 100s of Watts, 
where energy density is still an important consideration, but the 
scale is much larger than man-portable power generation where 
heat losses dominate [7]. Octane is used as a model fuel for gaso- 
line. This approximation is adequate from the point of view of 
mass balance and chemical equilibrium but would not be appro- 
priate for kinetic studies, and cannot capture catalyst poisoning 
effects. The above simplifications make the model independent 
of scale. 

The design specification is to produce a fixed amount 
of hydrogen from octane, water and atmospheric air (see 
Fig. 6) with maximal fuel energy density or maximal effi- 
ciency. Assuming an ideal fuel cell operating at T = 298 K and 
P = | atm a metric for the energy of hydrogen is approximately 
64 Wh (mol H2)~! [7]. Note that we use this upper bound on fuel 
cell performance, to isolate the conversion efficiency of the fuel 
processing step from the fuel cell performance. The fuel energy 
density is therefore defined as the equivalent energy of hydrogen 
produced, divided by the massflow of octane and water: 

fuel 64.N Hy out 


e — 


rav d 
z MW c,Hj3 NCsHig,in + MW4,0N420, in 


Since the hydrogen production Ny, out is fixed, maximizing 
the energy density is equivalent to minimizing the massflow 
of octane and water: 


max ete! e min(MWc,H is NCgHig + MWp, o Np,0). 


grav 


The energy conversion efficiency is defined as the equivalent 
energy of the hydrogen produced, divided by the consumed 
heating value of octane approximately 1251 Wh (mol CgHj3)~!: 


reac __ 64.N Hp out 
1251 NC;H;g,in 


The energy conversion efficiency is proportional to the selectiv- 
ity of hydrogen production, i.e., the molar flowrate of hydrogen 
produced divided by 9 times the molar flowrate of octane 
(NH, out) /(9NCgHj¢, in). Note that the selectivity can be greater 
than 1 from the conversion of water to hydrogen. Maximiz- 
ing efficiency is equivalent to minimizing the molar flowrate 
of octane. 


reac 


max n <> min(1251 NegHig,in) > min(NcgHijg, in): 


Q 
3 
y 
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Fig. 6. Flowsheet for the equilibrium-based case study. 


Note that for simplicity we do not account for the possibility of 
recycling, or further use of the excess water, carbon monoxide 
and methane. 

We assume a single-stage reforming reactor operating in the 
temperature range T = 800-1200 K and at ambient pressure 
which allows all species to be treated as ideal gases. We assume 
that at the reactor outlet the gases are at chemical equilibrium. 
The embedded chemical equilibrium leads to a bilevel program 
[29]. The outer program corresponds to the optimal operation, 
i.e., the choice of inlet composition that maximizes the effi- 
ciency or energy density. Constraints in the outer program are 
the energy balance and the specified hydrogen production. The 
chemical equilibrium is embedded as the inner program. Due 
to convexity of the inner program, this bilevel program can be 
reformulated as a single level nonconvex optimization problem 
[30]. For the derivation of the model used (see Appendix A). 

We solve the formulated optimization problem using BARON 
7.5 [31] available through GAMS 22.2 [32]. The key results of 
maximizing the energy efficiency versus maximizing the energy 
density are shown in Table 3. The two objectives lead to drasti- 
cally different operation as well as performance. Maximizing the 
energy conversion efficiency leads to a reduction of the energy 
density by approximately 60%. As expected, maximizing the 
energy conversion efficiency leads to a water rich feed, since 
the water molar flowrate does not affect the objective func- 
tion. The maximal energy density is achieved at a significantly 
higher temperature (1125 K) than the maximal energy efficiency 
(899 K). In the case of maximal energy density, the idealized 
partial oxidation reaction 


CsHıg + 402 > 8CO + 9H2 
is a good approximation for the computed overall reaction 
CgHig + 4.3502 

— 7.59CO + 8.56H2 + 0.0334CH4 + 0.366CO2. 


The computed overall reaction in the case of maximal energy 
conversion efficiency is not accurately described by one of the 
idealized reactions considered in [7] but rather from a combina- 
tion of these reactions. 

The above results were confirmed with Aspen Plus, version 
12.1 [33], using an RGibbs reactor with the optimal inlet compo- 
sitions and including the trace species that were neglected here. 
The discrepancy in the outlet composition is in the order of 1% 
while the discrepancy in the energy balance, due most likely to 
the constant heat capacity used, corresponds to an error in the 
temperature prediction of around 10 K. 

To demonstrate the tradeoff between the two objective 
functions we also calculate the Pareto-optimal curve. A Pareto- 
optimal or noninferior solution is a set of decision variable 
values, such that if we try to improve one objective, the other 
will be degraded [29]. Maximizing the energy efficiency and 
energy density are the two extremes of the Pareto-optimal set 
(assuming unique optimal solutions [34]). The use of Pareto 
optimization is well established in multiobjective optimization, 
e.g., [29]. The Pareto curve can be obtained via the ¢-constraint 
method, which in the case of two objectives is formulated as a 
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Fig. 7. Pareto optimal curve and inlet composition for the equilibrium-based octane reforming. 


Table 3 
Results for equilibrium-based case study 


Property Maximum Maximum energy 
energy density efficiency 

Octane relative molar flowrate 0.117 0.087 
Ncghj,in/ Nup, out 

Water relative molar flowrate 0. 1.18 
Nuy0,in/ Np; out 

Oxygen relative molar flowrate 0.508 0.462 
No, in/ NH, out 

Temperature ( K) 1125 899 

Energy density (Wh kg~!) 4800 2050 

Energy conversion efficiency (%) 37 50 

Hydrogen selectivity (%) 96 130 


single-parameter right-hand-side parametric optimization pro- 
gram [29]. We approximate this program by discretization with 
50 equidistant points in the energy efficiency space. The results 
are shown in Fig. 7. On the left side we show the Pareto curve, 
i.e., the maximal energy density as a function of the energy effi- 
ciency, and on the right side the corresponding relative molar 
flowrates of oxygen and water as a function of the relative molar 
flowrate of octane. The high consumption of octane corresponds 
to a high energy density and a low energy efficiency. 

The tradeoff between the two objectives is nearly linear. The 
increase of octane flowrate is associated with a slight increase 
in oxygen and a very pronounced decrease in water flowrate. 


4. Conclusions and future work 


For portable, and especially man-portable applications, the 
appropriate performance metric and design objective is energy 
density [6]. On the other hand, in large-scale power generation 
and energy conversion, the dominant objective is maximizing 
the energy conversion efficiency. We showed through analysis 
and case studies that these two objectives may lead to very dif- 
ferent design, operation and/or performance. In particular we 
showed that energy efficiency is not a suitable metric for select- 
ing, designing nor operating portable power generation devices, 
i.e., neither for choosing a given process among a list of candi- 
date alternatives nor for finding the optimal design and operation 
parameters of a given alternative. 


The models used largely neglect peripheral components such 
as valves and pumps, under the assumption that the influence 
of these components is not significant. Accounting for these 
components is of interest. We focused on fuel cell systems, 
but most observations are directly applicable to other systems 
transforming chemical to electrical energy, such as the use of 
thermophotovoltaic cells [35-37] or a microturbine driving a 
generator [38]. Also, studying hybrid electrochemical systems 
with a combination of a battery/super capacitor and a fuel cell 
is of interest. 
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Appendix A. Equilibrium-based model 


The chemical equilibrium of a set of species i € J contain- 
ing elements j € J is formulated as a nonlinear minimization 
problem with linear constraints and a strictly convex objective 
function [39] 

i out 


min X N; G(r) + RT In| — \+RT m N 
Nout me ome ' Pret oe Nk out 
s.t. 5 di iNi out = 5 @i Niis VYjJEJ 


ieJ ieJ 
Ni out = 0, Vie J 


where Nj in and Nj out for i € J denote the molar flow of species 
i into and out of the reactor, respectively. Note that the convex- 
ity of the objective function is with respect to the optimization 
variables (molar flowrates) and not with respect to parameters 
(temperature, pressure). The coefficients «;, ; correspond to the 
number of atoms of element j in species i. G;°(T) denotes the 
molar Gibbs free energy of the pure species i at the reference 
pressure Pref in the gas phase. Because of the logarithmic term, 
for implementation a nonzero lower bound needs to be imposed 
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on the molar flowrates N; out => € > 0 and this causes problems 
with trace species. We consider operation at ambient pressure, 
and take this as the reference pressure. 

To alleviate numerical difficulties, and based on the operat- 
ing conditions and previous calculations of chemical equilibria 
[7], we make some further simplifying assumptions. Nitrogen 
N2 is treated as an inert, neglecting the formation of nitrogen 
oxides and ammonia. This simplification is justified since the 
formation of nitrogen compounds is negligible from an mass 
and energy balance perspective and useful since it eliminates 
several species and an atom balance. Furthermore, complete 
conversion of oxygen and octane is assumed. This assumption is 
justified since, for the range of inlet compositions and tempera- 
tures considered, these species are only found in trace amounts 
at equilibrium. It is useful since the elimination of trace species 
makes the numerical behavior much more benign. Finally, we 
consider that the only hydrocarbon formed is methane. This 
simplification is again based on previous calculations [7], and 
is useful since it reduces the number of species. Since some 
species are eliminated, we distinguish between the set of inlet 
species Jin = {CgHig, H20, O2} and the set of outlet species 
Jout = {H20, H2, CO, CO2, CH4}. Note that nitrogen is elim- 
inated Nn, out = NN3,in = 3.79No,,in. The set of elements is 
given by J = {C, H, O}. 

It can be easily verified that all molar flowrates N; out are at 
nonzero value and therefore the inner program is equivalent to 
the following stationarity conditions (in addition to the element 
balances) 


min f(Nin) 
Nout, Nin; Tip i 


Assuming a constant heat capacity Cp and since the outlet is 
gaseous, 


G° (T) = HF (Tres) + Cpi(T — Tef) 


z HP (Trt) — Gi? (Tret) 
Tref 


+Cnin(=) 
In : 
s Tref 


where HË (Tief) is the ideal gas enthalpy of species i at reference 
temperature. 

Neglecting any heat losses, assuming that the outlet gases are 
at the reactor temperature, and that the inlet is at reference con- 
ditions (reference temperature Tref = Tamb = 298 K) the energy 
balance is given by 


XO Niou (HF (Tet) + Cri(T — Tret)) = X Niin Fi Tres). 
i€Jout i€Tin 
At reference conditions, water and octane are liquids at 
ambient conditions, while all other species are gases. 
Therefore Hi (Tret) = H? (Tret) for i¢{H20, CgHig} and 
Hi (Trt) = H? (Tret) — AH;"?(Tret) for ie {H20, CgHis}, 
where AH,” (Tef) is the vaporization enthalpy at the reference 
temperature. The enthalpies and Gibbs free energies in the gas 
phase and Tef are taken from [40], while the heat capacities 
and enthalpies of vaporization from [41], and are given in Table 
A.l. 

The optimal operation problem is formulated as the following 
single-level nonlinear program 


s.t.3.79Nop,in(C PN (T — Tret) + XO Niol HEt) + CPi — Tret) — X Niin HET) = 0 


i€ Jout i€ Jin 


XO aijNion— >> ajNiin =0, Yjeg 


i€ Jout i€Jin 


HP (Toet) — Gi? (Tref) 


H? (Tiet) + Cpi(T — Tret) r( 


N; 
+ RT In os 


Ferl (=) 
pi In | —— 
Tref ‘ Tref 


+ Qi jj = 0, Vi € Jot 
JEJ 


3.79Noy,in + kegan Nk.out << 


Np, out =1 

Niout = 9, Wi € Jout 
Niin = 0, Vi € Fin 
ujER, Vied, 


Ni, out 
NN3,out + De Tot Nk,out 


+ X oi, jj =0, Vi € Sout, 
jes 


G; (T)+ RT In 


where uj is the KKT multiplier associated with the balance of 
atom j. This formulation avoids the complementarity slackness 
conditions that are usually encountered when a bilevel program 
is transformed to a single-level program. 


where f(Nin) is one of the two equivalent objective functions, 
i.e., 


f(Nin) = NCsH;s,in 
for the energy conversion efficiency, or 
f(Nin) = MWeguig NCgHjg,in + MWe, 0 Nu, 0, in 


for the energy density. 
The formulated nonlinear program (NLP) contains noncon- 
vex functions and therefore global optimization methods are 
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Table A.1 
Physical properties 


Species HF (Tree) [J mol—!] AH? (Treg) [J mol~!] Gi°(Tret) [Jmol~!] Ch; Jmol! K7!] 
N2 0 N/A 0 30 
O2 0 N/A 0 32 
H2 0 N/A 0 29 
H20 —242 x 103 40 x 10° —229 x 10° 36 
co —110 x 10° N/A —137 x 10° 30 
CO; —394 x 103 N/A —395 x 103 47 
CH4 —75 x 10° N/A —51 x 103 50 
CsHig —208 x 103 42 x 103 —140 x 10° 150 


employed for its solution. It is relatively small with eleven vari- 
ables and two degrees of freedom and therefore can be solved 
easily to guaranteed global optimality with general purpose 
solvers. The degrees of freedom are two of the three inlet molar 
flowrates NcgHy¢,ins NH>0,in and No, in. The operating temper- 
ature is fixed by the energy balance. By the assumed constant 
heat capacity, the temperature could be eliminated from the set 
of variables by explicitly solving the energy balance, but this 
elimination is likely to make the optimization problem more 
difficult to solve. 
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